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Abstract 

Stationary vortex instabilities with wavelengths significantly larger than the thickness 
of the underlying three-dimensional boundary layer are studied with asymptotic methods. 
The long- wavelength Rayleigh modes are locally neutral and are aligned with the direc- 
tion of the local inviscid streamline. For a spanwise wave number P « 1, the spatial 
growth rate of these vortices is 0(/3 3 / 2 ). When P becomes 0(R the viscous correc- 

tion associated with a thin sublayer near the surface modifies the inviscid growth rate to 
the leading order. As p is further decreased through this regime, viscous effects assume 
greater significance and dominate the growth-rate behaviour. The spatial growth rate 
becomes comparable to the real part of the wave number when P = 0{R~ X ^ A ). At this 
stage, the disturbance structure becomes fully viscous-inviscid interactive and is described 
by the triple-deck theory. For even smaller values of P, the vortex modes become nearly 
neutral again and align themselves with the direction of the wall-shear stress. Thus, the 
study explains the progression of the crossflow- vortex structure from the inflectional upper 
branch mode to nearly neutral long-wavelength modes that are aligned with the wall-shear 
direction. 


1 INTRODUCTION 


Boundary-layer transition on a three-dimensional wing is often initiated by the so-called cross- 
flow instability that owes its origin to the inflectional nature of the boundary-layer profile along 
the crossflow direction (i.e., the direction that is orthogonal to the local inviscid streamline). The 
crossflow instability was first observed by Gregory and Walker (Gregory, Stuart, and Walker, 
1955) using napt.halene flow visualization, which indicated the vortex instability to be stationary 
relative to the airfoil surface. A theoretical explanation for this instability was given by Stuart 
in the same paper. He showed that the above disturbance modes are governed by Rayleigh’s 
equation, which is based on the boundary-layer profile along the direction of the disturbance 
wave-number vector. Stuart also went on to demonstrate that the Rayleigh problem permits 
unstable eigenmodes for each transverse wave number (i.e., wave number in the spanwise or 
the azimuthal direction, depending on the body geometry) which is smaller than some criti- 
cal value. The critical (i.e., neutral) wave number was shown by Stuart to correspond to that 
normal mode for which the inflection point of the directional mean profile coincides with the 
wall-normal location where the directional speed is equal to zero. 

The experiments of Gregory et al.(1955) established the rotating-disk boundary layer as a 
laboratory prototype for studying the crossflow-vortex instability. Viscous stability calculations 
for this flow were carried out by Malik (1986), who demonstrated the importance of retaining the 
Coriolis terms in the relevant (quasi-parallel) stability equations and also discovered the long- 
wavelength (or lower) branch of the neutral stability curve in the wave-number Reynolds-number 
plane. The asymptotic structure of these lower branch modes was simultaneously explained by 
Hall (1986). He showed that the lower branch vortices are aligned with the direction of the 
mean wall shear stress and, furthermore, are significantly influenced by the viscous and Coriolis 
forces. Hall also calculated the leading viscous correction to the inviscid (i.e., Rayleigh mode) 
neutral solution (generally referred to as the upper branch mode) which was investigated earlier 
by Stuart (Gregory, Stuart, and Walker, 1955). 

If the Rayleigh mode calculation for the rotating-disk boundary layer is continued from 
the upper branch wave number toward smaller values of the azimuthal wave number, then the 
radial amplification rate increases initially, reaches its maximum at some finite wave number 
and, finally, falls back to zero again as the azimuthal wave number tends to zero (Choudhari 
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and Streett, unpublished results, 1990). However, somewhat surprisingly, the limiting form of 
the long-wavelength modes corresponds to vortices that are aligned with the local free-stream 
direction, against the viscous Coriolis-balanced vortices of Hall (1986), which are aligned with 
the direction of the mean wall shear stress. A similar observation was recently made by Blackaby 
and Choudhari (1993) in the context of the Falkner-Skan-Cooke (FSC) family of boundary layers, 
in which the Coriolis forces are absent altogether. Because Coriolis forces are also absent from 
the boundary-layer flows encountered in typical aeronautical applications, an investigation of 
the properties of the long-wavelength Rayleigh modes is deemed important; this investigation is 
the purpose of the present article. 

In section 2, we show that the long-wavelength inviscid vortices are neutral to the leading 
order and that their relatively small amplification rate is determined by the local sweep angle 
and the momentum integral associated with the jetlike crossflow profile. The latter dependence 
is characterized analytically, and the similarity of the result to that obtained earlier by Drazin 
and Howard (1962) in the case of unbounded, two-dimensional jets and wakes is pointed out. 
The stabilizing effects of convex surface curvature on the long-wavelength crossflow instability 
are considered in section 3, wherein we elucidate the asymptotic behaviour of the relevant 
neutral branch found by Blackaby and Choudhari (1993). Perhaps the most significant difference 
between the case of free-shear flows examined by Drazin and Howard (1962) and the boundary- 
layer case examined herein is related to the earlier onset of viscous effects in the present problem. 
The role of viscosity in modifying inviscid growth-rate behaviour as the transverse disturbance 
wave number /? decreases through the successive regimes of /? = 0(R 1|/7 ), = 0(R ^ 4 ), and 

/? = 0(R ~ 3 / 7 ) is discussed in sections 4.1, 4.2, and 4.3, respectively. Concluding remarks based 
on this study are presented in section 5. 

2 ANALYSIS 

For a three-dimensional mean boundary layer, the profiles along a pair of suitably chosen Carte- 
sian axes x* and z* along the surface are assumed to be given by V* cos 0 SW u{Y\ x, z) and 
V* sin 9 SW w(Y] x, z), respectively. (See figure 1.) Here, x = x*/(* and z = z* / 1* are the nondi- 
mensional coordinates based on the length scale £* over which the boundary-layer flow evolves; 
Y = y* / L* is the wall-normal coordinate scaled by a typical thickness L* of the boundary layer; 
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and 6 SW is the local sweep angle, so that u(V') and w(Y) represent the normalized profiles at 
station (x,z) (i.e., limy _ 00 n(y) == limy^oo w(Y) = 1). Without any loss of generality, we will 
assume that — 7t/2 < 9 SW < 7r/2, cos9 sw — 0(sin 0 SW ) (so that the local free stream is not aligned 
with either the x or the z axis) and that u ^ w for all Y (so that the base flow is genuinely three 
dimensional). In the case of swept- wing boundary layers in low-speed flows, it is usually true 
that the velocity profile along the direction of the inviscid streamline (i.e., cos 2 9 SW it + sin 2 9 SW w) 
has a monotonic behavior in Y. Also, the crossflow profile (= cos 9 SW sin 9 SU) (u — w )) reaches its 
maximum magnitude within the boundary layer and exhibits a single inflection point between 
the location of the above maximum and the edge of the boundary layer where the crossflow 
magnitude tends to zero. In order any possible complications associated with more intricate 
boundary-layer profiles, we will restrict our attention to profiles that satisfy the above proper- 
ties. 

As stated in the introduction, the objective of our analysis is to trace the evolution of the 
crossflow- vortex characteristics as the wavelength of the vortex is increased gradually, starting 
from its value for the inflectional, upper-branch mode. It is shown in sections 2 and 4.1 below 
that in the first two asymptotic subregimes encountered during this process, the disturbance 
structure involves an asymptotic subregion near the edge of the mean boundary layer, where 
the crossflow magnitude is asymptotically small. In order to simplify the analysis, we assume 
that the mean-flow vorticity decays at an exponential rate as Y —> oo. That is, 

U ~ 1 - Cu exp [~k u {Y - Ys)] (la) 

and 

w ~ l - C& exp [-k w (Y - Y 6 )] (16) 

where (7^,(7^, /c u , k w , and Ys are known constants. In order to ensure that the variations in 
both u and w are significant in the abovementioned edge region (where Y >> 1) we assume 
that k u = k w = k. However, we shall demonstrate that the specific form of the decay in mean- 
flow vorticity is unimportant in determining the leading-order growth rate of the crossflow 
eigenmodes. Hence, the results obtained herein are equally applicable to boundary-layer profiles 
that have, for instance, a faster-than-exponential decay of vorticity into the free stream (e.g., the 
self-semilar Falkner-Skan-Cooke boundary layers). To minimize the influence of any extraneous 
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parameters, the curvature of the surface within the region of interest is also assumed to be 
negligible to begin with. 


2.1 Inviscid Treatment 

Now consider the linear, normal-mode disturbance with a spanwise wave number equal to ft* = 
ft/ L*,ft (> 0) = 0(1). As shown by Hall (1986), the dynamics of this disturbance is governed 
by the Rayleigh’s equation up to 0(/? _1/3 ), where R (» 1) is the Reynolds number based on 
the velocity V* and length L* . Thus, the perturbation quantities in the boundary- layer region 
expand as 

= exp [i(aX + ftZ)], (2) 

where the parameter A(<< 1) characterizes the infinitesimal disturbance amplitude; X — x */ L * 
and z = z* j L* are the “fast” coordinates scaled by the boundary-layer thickness L*\ and O' is the 
complex wave number along the X direction, which is the eigenvalue to be determined during 
the course of the analysis. The profiles [tt(V'), v(Y), w(V'), p(K)] represent the corresponding 
eigenfunctions that illustrate the variation of the respective disturbance quantities along the 
wall-normal direction. In particular, u(V) is given by the solution to the Rayleigh’s equation 

Crv = (a cos 9 SW u + ft sin 9 SW w) jv" — (« 2 + /3 2 )u] — (o' cos 9 SW u" + ft sin 9 SW w”)v = 0, (3a) 


which also satisfies the homogeneous boundary conditions 

v(0) = 0 and v{Y ) — ► 0 as Y — > oo 


(36, c) 


The primes in (3a) denote derivatives with respect to the wall-normal coordinate Y ; unless 
specified otherwise, the same convention will be valid throughout this paper. 

The objective of our long-wavelength analysis is to determine the asymptotic behaviour of 
the solution to the eigenvalue problem (3a-c.) when the transverse wave number ft tends to zero 
(i.e., R~ T « ft « 1 for any r > 0). As mentioned in the introduction, numerical calculations 
suggest that a ~ —ft tan 9 SW as ft — > 0 because the vortices tend to align themselves with the 
local direction of the inviscid streamline. As a result, only two disparate length scales exist 
in the long-wavelength problem, namely the boundary-layer thickness L* and the disturbance 
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wavelength 2ir / 0* . Therefore, one would intuitively expect the disturbance structure to split into 
at least two distinct subregions along the surface-normal direction, namely those that correspond 
to Y = 0(1) and Y u = 0Y = 0(1), respectively. However, the form of the asymptotic series 
relevant to either of these subregions cannot be predicted on an a priori basis because of the 
singular nature of the limit (3 — > 0. 

The correct form of the asymptotic expansion can be determined with the benefit of some 
hindsight based on a knowledge of the properties of Rayleigh’s equation, although in our deriva- 
tion we were also guided by the closely related analysis of Papageorgiou and Smith (1989) on 
the long-wavelength stability of two-dimensional wake profiles. In view of this, the issue of jus- 
tification for the form of asymptotic expansion used is postponed until later; the reader is asked 
to assume until then that the expansion for both the eigenvalue a and the eigenfunction v(Y ) 
across all the relevant subregions proceeds in integral powers of (3 (at least up to the order 
required to determine the leading-order growth rate). Given this form of the expansion, we 
find that in addition to the two subregions already mentioned (i.e., Y = 0(1) and Y u = 0(1)) 
a third, intervening region that has the same thickness as that of the boundary layer and is 
located near the outer edge must be considered at a distance of 0(log 0) 1 ^ 2 from the surface, 
where the magnitude of the crossflow velocity is O(0 1 ^ 2 ). 

We now begin the solution procedure by rewriting the Rayleigh’s equation in a form that is 
suitable for long-wavelength analysis. Thus, (3a) is expressed in the form 

Crv = (Cw + 0 2 M) v = 0. (4 a) 

The two operators £yv and M. are defined by 

C„v s [w 2 (^)'J' (46) 

and 

M.v = — (1 + tan 2 6 sw a 2 ) W v, (4c) 

wherein we have introduced new notation in the form of a scaled eigenvalue parameter 

at = a/0 tan 9 SW [Ad) 

and the normalized directional profile 

W(K) = m/(F) + w(Y) . (4e) 
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Ill anticipation of the result quoted earlier, a expands in the form 


a ~ — 1 -f /5 1//2 «i + fiot 2 + 3 + f3 2 ot4 + ••• (/? — ^ 0) . (4/) 

Therefore, the operator norms for both £yy and M continue to be of 0(1) in the main boundary- 
layer region during the above limiting process. 

For Y = 0(1), (4f) implies that 

W(Y) ~ Wo + u + (3 a 2 u + /? 3/ ' 2 a 3 u + (3 2 a 4 u + ... [Wo(K) = w(Y) - u(Y )] . (5a) 

Accordingly, the eigenfunction u(V") can be expanded as 

v(Y) = pv 0 (Y) + /3 3/2 v,(Y) + (3 2 v 2 (Y) + p 5/2 v 3 (Y) + /3 3 v 4 {Y) + ... , (56) 

where the leading order perturbation has been taken to be 0(/9) so that the magnitudes of other 
disturbance quantities u y w , and p remain 0(1) even in the limit /? — ► 0. If we substitute the 


expansion (5b) into (4a) and use (4f) and (5a) with j3 — ► 0, we find that u n ,n — 1 — 4 satisfy 

Cw 0 v o = 0, v o (0) = 0; (6a, b) 

Cw 0 vi = -a A CuV 0 , vi(0) = 0; (7a, 6) 

£wo u ‘2 = ^~ a 2 ^uVo — o^\CuV\^ ^(0) = 0; (8a, b ) 

Cw 0 v 3 = -ot 3 CuV 0 - ot 2 CuV 1 - ot\CuV 2 , v 3 (0) = 0; (9a, b ) 

and 

£yV 0 V 4 = “^OL^CuVq — Ot 3 CuV\ ~ OL 2 C'u v 2 ~ — v 4(0) = Oj (10a, b) 

where 

C Wq v = VV 0 2 (p^-) , CuV = u 2 (^) , and M 0 = -sec 2 0 3W VV 0 . (10c, d,e) 


The behaviour of v n (T), n = 1—4, when Y — ♦ oo has not been specified at this stage, because 
it is determined by asymptotic matching with the solution from the adjacent outer region. 

The general solution to the homogeneous, long-wavelength Rayleigh equation (6a) is given 

by 

u 0 (Y) = C a v a (Y) + C b v b (Y), (11a) 
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where 

v a (Y) = W 0 (Y) ~ Wo'(O) Y + 0{Y 2 ) ( 116 ) 

and 

^= w °£w~vm +0{Y ' ogY) - (llc) 

Observe that 

^(n-Wo'^r + ocr 2 ), (ii d) 

and 

” s(5/) ~vvbW + 0(l,|ogr) < lle) 

as Y —* 0, and 

v A (Y) ~ Qwe~« Y , (Q w = (C Q - C w )e KY *) , (11/) 

and 

MY) ~ ^ e_ “ r ’ < u ») 

as y — ► oo. In view of (lid) and (lie), we must set Cb = 0 to satisfy the surface boundary 

condition (6a); the leading order eigenfunction is, therefore, given by 

= WoOO, (12) 

wherein Ca — 1 without any loss of generality. Note from (Ilf) that the above solution already 
satisfies the outer boundary condition (3c) on v because a = 1 has been set to the leading order. 
The higher order solutions v ni 7 1 =1—4 are now given by 

V\ — a\ u, V 2 — ot‘ 2 u, (13a, 6, c) 


and 

V 4 = <^ 4 u + sec 2 9 SU) Wo f — 2 / Wo^V^dl'i, (13d) 

Jo Wo Jo 

respectively, and the coefficient of va has been set equal to zero in each case to satisfy some 
(arbitrarily chosen) normalization constraint. Because u — > 1 as Y — > oo, none of the solu- 
tions from (13a) through (13d) satisfy the original outer boundary condition (3c). Although 
Ui through U 3 at least retain their 0(1) magnitude as Y — ► oo, the second term in v 4 becomes 
exponentially large there. For large values of F, (13d) can be rewritten in the form 

v 4 - a 4 u + Io(O sw ,Wo)v b {Y) + sec 2 0 SW VV 0 / r- 7 ^ f W 0 2 dY 2 dY l , (14a) 

Jo Wo Jo o 
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where the functional I 0 in the second term is defined as 

fOO 

Io(0 sw , Wo) = sec 2 e sw / Wo 2 dY (14 b) 

J o 

so that the third term becomes exponentially small as Y — ► oo. 

The ordering within the asymptotic series (5b) breaks down when /? 1/,2 iq becomes comparable 
in magnitude to the leading-order term v Q (i.e., when Wo(K) = 0( /3 1 / 2 )). In view of the assumed 
behaviour (la, b), this breakdown occurs when Y = O [— T-log(/?)] (i.e., in the region close to 
the outer edge of the mean boundary layer). (See figure 1.) The new coordinate is introduced 
as 

Y = y + iog(/J) = 0(l), (15) 

2k 

which is appropriate for this region. The vertical velocity perturbation v is expanded in the 
form 

v = /? 3/2 5,(T) + pv 2 (Y) + p 5/2 v 3 (Y) + ... , (16) 

which is suggested by the outer behaviour of the main-deck solutions (12) and (13a-d). If we 
note that the coefficient function W expands as 

W ~ /? 1/2 Wi + pm + /? 3/2 W 3 + ... (17a) 


in the edge-layer region where 



VVi = Qyy e ^ “I" {Qw — {Cu “ 

(176) 

and 


W„ = QuCtn-ie K? + On for n > 1 ( Qu = C ii e KYi ') , 

(17c) 

then 

we find that the perturbation quantities v n , n — 1 - 3 are governed by the equations 


o 

II 

(18a) 


CyjV* — QL\C,'fy 2 V\i 

(184) 

and 


^VVi ^3 "™ — a '2^VV 3 — A* A 

(18c) 


The operators £ V y ] , £yy 2 , and C are defined analogously to those in (10c, d). 
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Note that the expansion (16) has been truncated at the third term (i.e., at 0(fi 5 / 2 )) because 
the latter will match the leading contributions from both v 3 and v 4 from the main boundary- 
layer region. The solutions for the edge-layer perturbations, for which the inner limit as y — * oo 
matches the outer limit of the main-deck solutions (12) and (l.‘3a-d), are given by 

= v a (Y), (19a) 


v -2 = VV 2 (F), (196) 

and 

n 3 = W 3 (Y) + C 3 v a (Y) + I 0 vb(Y), (19c) 

where 

v a (Y) = W x {Y) (190 


and 


vb 


(?) = vv, / 

Jo 


y dY 

vvf 


(19c) 


analogous to (11b, c); C 3 is a constant that does not need to be determined for this analysis. 
Because the solutions (19a) through (19c) can also be obtained by expressing the main-deck 
solution in terms the edge-layer variable Y and reexpanding in powers of /?, the edge layer plays 
a passive role in determining the overall dynamics of the long-wavelength disturbance. 

Because YV\ —► a\ as Y — > +oo, neither Ct, n 2 , nor v 3 satisfy the original outer boundary 
condition (3c). Moreover, v 3 exhibits a linear growth 


- /o vy 

y 3 > ag 

^ 1 


OO, 


( 20 ) 


which makes the asymptotic series (16) invalid at sufficiently large values of Y. This linear 
growth of the eigenfunction magnitude is cut off within the potential-flow region (or the upper 
deck), where Y u = (3Y = 0(1). As we will observe below, the leading-order eigenfunction within 
the potential-flow region has a simple, exponentially decaying character. The zeroth-order (i.e., 
the constant) term from the small- V u Taylor expansion of this eigenfunction matches the outer 
behaviour of v x from the edge-layer region. The coefficient of the linear term from the above 
Taylor expansion now matches with the slope of the linearly growing component (20) from the 
edge-layer solution to provide the dispersion relationship for the eigenvalue correction a,, which 
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determines the spatial growth rate of the long-wavelength crossflow vortex. As seen from (20), 
the slope of the outer limit of v 3 is completely determined by the properties of the mean-flow 
profiles within the main boundary-layer region and does not depend on the assumed form of 
vorticity decay as Y — ► oo. Thus, the dispersion relationship obtained will actually be valid for 

a much wider class of boundary-layer profiles. 

Only the leading-order solution must be considered within the upper deck region, wherein 



W ~ /3 ,/2 « i + - 

(21a) 

and 


V ~ + .... 

cm 

The Rayleigh’s equation (4a) 

now reduces to 



d 2 v i 2 

M/9 sec 0 8W v i = 0 
a i 1 

u 

(22) 


to the leading order; the solution to (22), which satisfies the outer boundary condition (3c), is 
given by 

vi =u 1 (0)e- sec<? '*"’ r “. (23) 

A match with the leading-order solution from the edge layer requires that 

«i(0) = a, (‘24) 


and the linear component of the Taylor expansion of V\ about V u — 0 matches onto (20) provided 
that 

(25) 


h a 

— = — sec0 su ,c*i, 

«i 


where I 0 {6 sw ,u,w) is defined by (14b). Equation (25) yields the dispersion relation 

Io \ ,/a 


“ ,=±, (sdc) • 


(26) 


If we restore the above result to the streamwise wave number o, then we have the asymptotic 
expansion 


« ~ -p tan fl„ ± ip 3 ' 2 1 /“ (<o - u)‘ dvX" 

COS 3 / 2 v sw Uo J 


as /3 —*■ 0. 


(27) 
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Equation (27) shows that the first correction to the leading-order eigenvalue is purely imag- 
inary. Accordingly, it corresponds to a relatively small rate of amplification or decay along the 
X direction, depending on a. the sign of sin 9 SW and b. whether the minus or the plus sign is 
chosen in (27). For sin0 jur > 0, the choice of minus sign leads to a spatially unstable normal 
mode, whereas the plus sign needs to be chosen if sin0 sw < 0. The appearance of two normal 
modes (both of which have the same wave number in the X direction) is not surprising in view of 
the asymptotic expansion (4f), which proceeds in integral powers of /7 1 / 2 . Equation (4f) implies 
that, even though d^jda — * 0 as /? — » 0, d 2 /3/da 2 becomes singular in this limit. This behaviour 
explains the observed coalescence of the unstable and stable modes in the long- wavelength limit. 
The results of figure 2(a) will confirm that the unstable mode does indeed describe the long- 
wavelength behavior of the unstable Rayleigh solution that is obtained numerically, in general. 
It should be noted that the unstable modes is physically relevant only if it can be realized as the 
time-asymptotic limit of the corresponding initial boundary value problem. (In other words, if 
the underlying boundary-layer flow is convective unstable near the location of stability analysis 
(Huerre and Monkewitz, 1990).) We will assume that the normal mode is physically realistic 
throughout this paper. See the remarks by Choudhari (1994) in this context. 

As is common with long- wavelength instabilities, the amplification rate implied by (27) is 
independent of the details of the velocity profiles that correspond to the underlying boundary- 
layer flow and is completely determined by both the local sweep angle (measured relative to 
the direction of amplification) and the momentum integral associated with the jetlike crossflow 
profile. A similar result was shown to apply in the case of two-dimensional wakes and jets 
in the classic paper by Drazin and Howard (1962). They provided a physical explanation for 
the result with an estimate of the centrifugal force associated with the curved streamline that 
results from particle displacement induced by the disturbance. The paper by Papageorgiou and 
Smith (1989) placed Drazin and Howard’s analysis in a proper asymptotic framework, and we 
have shown an analogous structure to be applicable to the stationary crossflow instability of 
three-dimensional boundary-layer flows. During this process, we have also indicated why the 
result (27) is independent of the specific form of vorticity decay that is associated with the 
underlying boundary layer. It is also interesting to note that in spite of the spatial growth rate 
being asymptotically small with respect to the real part of the wavenumber, the long- wavelength 
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inviscid mode does not have a critical layer inside the boundary layer. 

Now that a self-consistent solution has been derived for the long-wavelength Rayleigh prob- 
lem, we can justify why the eigenvalue expansion (4f) was chosen to proceed in powers of /3 1 / 2 
rather than, for example, in powers of 0 2 as a regular perturbation would have indicated. There- 
fore, suppose that the eigenvalue expansion (4f) proceeds in integral powers of 0 l ' J . That is, 

OO 

a ~ — 1 + (28) 

j = o 

where J is an integer to be determined during the course of the analysis. This section demon- 
strates that the leading term of the corresponding expansion for the eigenfunction is a simple 
displacement of the jetlike crossflow profile in the wall-normal direction. As indicated by (7b), 
the first-order correction to this displacement behaviour is caused by the O(0 l/J ) deviation in 
the directional profile W from its jetlike behaviour Wo, which is forced by the eigenvalue correc- 
tion c*i in (28). The discussion that follows (19e) implies that the above correction term can be 
continued across the edge layer without any change in its order of magnitude and that it rapidly 
tends to an O(0' l,J ) constant near the outer part of the edge layer. On the other hand, the O(0 2 ) 
correction term, which in part is forced by the operator M in the Rayleigh’s equation, gives 
rise to the singular component (i.e., vb) in the eigenfunction. This part exhibits an exponential 
growth in the main deck, which leads to a magnification by 0(0 1//2 ) in its order of magnitude 
as Y — > oo. The exponential growth is cut off within the edge layer and is converted to a linear 
growth that amounts to an additional magnification by 0(0 J ) in the magnitude of v& as the 
wall normal location varies across the edge layer. The 0(/3 2— ^ — 1 ) linear component produced in 
this manner is commensurate with the 0(0$) constant behaviour at the inner edge of the upper 
deck only if 

(2-i)- 1 ^ (29) 

(i.e., if J = 2 as assumed to this point). 

To verify the predictions of the long-wavelength asymptotic theory and to estimate the range 
of wave numbers in which the relevant approximation is reasonable in a numerical sense, the 
asymptotic predictions based on (27) must be compared with those obtained from a numerical 
solution of the Rayleigh’s equation. Of course, it should be borne in mind that the primary 
purpose of asymptotics is not to provide numerical results, but to pinpoint the dominant physical 
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processes that determine the structure and the growth properties of the unstable disturbance. 
For convenience, we have opted to illustrate the above comparison in the case of a self-similar 
FSC profile with the Hartree parameter equal to 0.6 and a local sweep angle of 9 SW = 45°. 
To enable the reader to form an overall perspective, the results of the Rayleigh calculation 
throughout the region of instability are shown in figure 2(a). The small wave-number portion of 
this graph has been replotted in figure 2(b); the additional dashed curve indicates the asymptotic 
result. Note that in view of the nearly linear behaviour of the Real(a) curve in figure 2(a), the 
corresponding long-wavelength prediction remains accurate throughout the region of crossflow 
instability. This behaviour is a result of the narrow directionality of the stationary crossflow 
vortices. On the other hand, the graph for Imag(a) in figure 2(b) indicates that the spatial 
amplification rate begins to deviate from its square-root behaviour at small values of the wave- 
number parameter /?. 

For the self-similar family of FSC boundary layers, the profiles u and w are dependent 
only on the Hartree parameter /?#; for flu = 0.6, the normalized crossflow-momentum integral 
f£°(w — u) 2 dY has a value of 0.039. Note that for the self-similar FSC profile, Y corresponds to 
the similarity coordinate y*[V*cos 6 sw /{2 — (3 h)v*x*] x I 2 , i.e., L* = [V*cos 0 sw /(2 — /?//)i/*x*] -1 / 2 . 
For future reference, figure 2(c) shows the variation in the latter quantity as a function of the 
Hartree parameter i3h- 

A similar comparison for the real and imaginary parts of a as functions of the sweep angle 
0 SW at a fixed wave number of (3 = 0.01 is indicated in figure 2(d). Because (3 = 0.01 represents 
a small value for the wave number parameter, the asymptotic prediction for Real(a) is virtually 
indistinguishable from the corresponding numerical result based on the Rayleigh calculation. 
The graph for Imag(a) shows that the long-wavelength prediction for the spatial growth rate is 
satisfactory throughout the range of sweep angles considered; however, a slight loss of accuracy 
is observed when 9 SW —y 90°. 

3 Effect of Convex Streamwise Curvature 

For the analysis in section 2, we had assumed that the curvature of the body surface was negligi- 
ble. However, cases of practical interest (Collier and Malik, 1987) exist for which the curvature 
is strong enough to modify the growth rate of the stationary Rayleigh modes. The effects of sur- 
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face curvature are likely to be most significant (in a relative sense) when the associated crossflow 
instability is weak (i.e., when 0 SW — > 0). This limit corresponds to a mean flow that is predomi- 
nantly along the positive X direction across the entire boundary layer. However, note that the 
existence of crossflow instability depends on a balance between the streamwise and spanwise 
convection terms in the governing equations. This balance, which is reflected through the com- 
bination a cos 9 sw u + >3 sin 6 sw w in ( 3 a), implies that when 0 SW -+ 0 , a — ► 0 in such a way 
that a remains finite. A further balance between the convective and the centrifugal terms in the 
normal-momentum equation (not shown in this paper) indicates that the normal-mode growth 
begins to be comparably affected by both crossflow and surface curvature terms when the mag- 
nitude of the local sweep angle 0 SW has fallen to an 0{R~ 1 ^ 2 k x ^ 2 ) level. Here, k = L* fr (<< 1) 
denotes the nondimensional curvature of the surface based on Tq, which is the local radius of 
curvature along the X direction; we have assumed the other Gaussian curvature of the surface 
(termed the “transverse” curvature in literature) to be zero for convenience. 

The parameter region that involves the combined effects of crossflow and surface curvature 
may also be approached from the other extreme in which the starting point of the analysis is the 
Gortler instability of a two-dimensional boundary layer, and the crossflow parameter (i.e., the 
local sweep angle) is allowed to increase from its initial zero magnitude with the transverse wave 
number f) fixed at some 0(1) value. This approach was taken by Bassom and Hall (1991a), who 
extended the pioneering studies of Hall (1982, 1985); Denier, Hall, and Seddougui (1990); and 
Timoshin (1990) to include three-dimensional boundary layers. As shown by Bassom and Hall, 
the inviscid instability under the combined effects of weak crossflow and surface curvature is 
governed by the Rayleigh’s equation for streamwise-elongated vortex structures. This equation 
contains one additional term that characterizes the influence of centrifugal forces associated with 
the surface curvature. After the introduction of the scaled sweep angle 8q = ^72 where the 
Gortler number G is defined as G = 2 Rk (>> 1), equation ( 2 . 8 a) of Bassom and Hall (1991a) 
can be expressed in the form 

X + » = 0 . (30a) 

Recall that the operators and M in (30a) were defined earlier in (4b) and ( 4 c) and the 
parameter \ equal to ±1 depending on whether the curvature is concave or convex, respectively. 

Because the primary interest of Bassom and Hall’s investigation was related to Gortler 
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instability, they restricted their calculations to the case of concave surface curvature. Their 
results confirm the findings of Malik and Poll (1985) and Collier and Malik (1988), which show 
that this type of curvature has a destabilizing effect overall. Blackaby and Choudhari (1993) 
focused their interest on the neutral solutions to the Bassom and Hall equation. They considered 
both concave and convex curvatures and established a connection between the corresponding 
neutral curves by drawing an analogy to the problem of the stability of stratified shear flows. 
One of their findings, which is of interest to us in this paper, is the existence of a neutral curve 
that involves a long-wavelength branch in the case of a convex surface. This study showed that 
the weakly three-dimensional flow over a convex surface is completely stable when the crossflow 
magnitude is below a certain critical value. As the crossflow parameter exceeds this critical value, 
disturbances in a finite band of wave numbers become unstable. The bandwidth of unstable 
disturbances increases with the amount of crossflow applied. In the limit of large crossflow 7 
magnitudes, the right branch of the neutral curve (which bounds the region of instability in 
the 0 — 0 SW plane) approaches the inflectional neutral mode of the Rayleigh solution (Gregory, 
Stuart, and Walker, 1955); the left branch tends to the long- wavelength limit (0 — ► 0) that was 
considered in section 2. 

We now illustrate a simple modification to the analysis in the previous section to account 
for the stabilizing effect of convex curvature. This modification yields a simple result for the 
long- wavelength asymptote of the neutral curve found by Blackaby and Choudhari (1993). As 
indicated in section 2, the amplification rate of long-wavelength crossflow 7 vortices is determined 
by the O(0 2 ) term in equation (4a), which involves the operator A4. Clearly, in order for the 
curvature term to nullify this effect, we must have 0 O ~ ^ (0 O o = 0(1)) as 0 — » 0. In that case, 
the asymptotic expansions for the different subregions involved, as w 7 ell as the relevant solution 
procedures, remain completely analogous to those described in section 2 above. The first major 
difference occurs in solving for V 4 in the main-deck region, which now satisfies 

Cw q V 4 = -a 4 C u v 0 - a 3 jC u vi - a 2 CuV 2 - a^C^v 3 - Mgo v o v 4 (0) = 0 (31a) 

instead of (10a, b), where 

Mgo v o = ~ ^VVq 2 + ^2 ^ • (31^) 
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Equations (31a, b) lead to 

,y 1 r Y i / roo 1 \ 

u 4 = a 4 u + Ici(Qoo, Wo)vs(E) 4- W 0 / / MaoWodY^dYi f Ig - / W 0 2 dy' - 75 - ) , 

-'O Wo J oo \ Jo ^00/ 

(32) 

which is analogous to (14a). The solutions in the edge-layer region and the upper deck remain 
the same, except for the replacement of Iq by Iq. This result leads to the dispersion relation 

c*i = ±i Ia 1/2 (33) 

that determines the amplification rate of the vortex. The neutral curve, which corresponds to 
the imaginary part of equal to zero, is then given by 

Io = 0 . 


That is, 


do 


P { 2 / 0 °° Wo 2 dY)' /2 


as /3 — > 0 . 


(34) 


Note that the large-crossflow behaviour implied by (34) is more relevant in laminar flow control 
than the 6 0 = 0 ( 1 ) case, in view of the strong three dimensionality that is involved in such 
applications. A comparison of the analytical result above with the numerically obtained neutral 
curve by Blackaby and Choudhari (1993) is shown in figure 3. The mean flow chosen for this 
comparison corresponds to the FSC profile with a Hartree parameter of /?// = 1 / 2 . (See the paper 
by Blackaby and Choudhari for details.) Figure 3 shows that the long- wavelength asymptote 
closely approximates the behaviour of nearly the entire left-hand branch of the neutral stability 
curve in the 9 0 — plane. 


4 Viscous Effects 

4.1 Viscous-Modified Inviscid Regime 

Now consider the effect of viscosity on the long-wavelength solution derived in the previous 
subsection. Of particular interest is whether or not the viscous effects can nullify the small 
inviscid growth rate of the long-wavelength vortices to give rise to a lower branch of the neutral 
stability curve. The answer to this question has important consequences because transition 
prediction methods of today utilize the viscous, quasi-parallel framework to determine the local 
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amplification rate, which is then integrated from the lower to the upper branch station in 
order to estimate the location of transition onset (Bushnel! and Malik, 1987). In addition, the 
lower branch properties are particularly important in predicting the receptivity 1 of the three- 
dimensional boundary layer. 

Viscous action within the inviscid regions of Y = 0(1) and larger can only modify the vertical 
velocity perturbation found in section 2.1 through a correction term of 0[(/3/?)“ 1 ]. However, 
a much larger correction (of 0[((3R)~ 1 ^ 3 }) is provided by the viscous sublayer close to the wall, 
which is required so that the disturbance satisfies the no-slip condition at the surface. The X 
and Z components of the inviscid slip velocity associated with the main-deck solution (12) are 
given by 


?z(0) = 


i-iM. 

tan 6sw 


+ 0(jl' r ‘) 


(35a) 


and 


ui(0)=iw'(0) + O(/3 ,/2 ) . 


(35 b) 


The slip velocity given by (35a, b) is brought to zero within the viscous sublayer (or the 
lower deck), where Y] = [fiR^^Y = 0(1). The disturbance quantities within this region are 
expanded in the form 


(u,v,w, P )= + ... . 


( 36 ) 


By substituting this expansion into the continuity and momentum equations and allowing (3R 
oo with Yi — 0(1), we are led to the linearized three-dimensional boundary-layer equations 



i(W 0 tan 9 sw Uo) + ... =0, 

a / / 

(37a) 

d 2 U 

:W 0 (O)cos 9 sw Yi ] U 0 + cos 9 sw u'( 0) V 0 = -i tan 9 SW P 0 + 2 , 

a Y\ 

(376) 


dPn 

(37c) 


O 

II 

5 ^ 

[?W(/(0)cos 

9 S wYi] Wo + sin 9 sw w (0) V 0 = iP 0 + 2 . 

d Yi 

(37rf) 


*The term receptivity refers to the various physical mechanisms through which the external disturbances 
generate the unstable eigenmodes of a boundary-layer flow (Morkovin, 1969). 
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The boundary conditions for the lower deck velocities correspond to 


U 0 = V 0 = W 0 = 0 at Y, = 0 


(38a, b, c ) 


and 

Uo — u(0), Wo = u>(0) as V/ — oo (38d, e) 

to match with the inner limits (35a, b) of the main-deck solutions. 

If the equations (37a-c) are solved along with the boundary conditions (38a-e), then the 
vertical-velocity perturbation Vo is given by 


Vb = 


Vo, 


C, £ Ai(i)di - Ai'(Q + Ai'(Q) 


(39 a) 


Ai'(0) |/ 

where Ai and Ai' denote the Airy function of the first kind (Abramowitz and Stegun, 1964) and 
its derivative, respectively. We have set 


3Az (0)Wo (0) r™> ><(\\ ■ a Al*v 

Wo = 1 777^ Hffi ’ C = [^o (0) sin 6 SW ] Y t 

[iW 0 (0) sin 9 SW ] ’’ 


(39&, c) 


where the principal branch has been chosen for the cubic root in the definition of Wx> and £ so 
that as Y[ — » oo 

Vo ~ W o '(0)V/ + Wo- (40) 

The constant term within the solution (40) for the vertical velocity perturbation gives rise to a 
correction to the main-deck solution (12)-(13d) that is of the form 


^vis 


= - 0 ({ 3 R)-V 3 WonV(o)Mr). 


(41) 


This viscous correction first affects the spatial growth rate of the long-wavelength vortex when 
it becomes comparable to the [Pv 4 inviscid term, which (shown in section 2) determines the 
inviscid amplification rate by producing a similar vb component as Y — » 00 . The viscous and 
inviscid contributions to the coefficient of the vb component are of the same order of magnitude 
when 

/ 3 — R~ A / 7 j3 0 (/?o = 0(1) as R — > oo), (42a) 

in which case the dispersion relationship (27) is altered to 

a ~ -/T 1/7 A) tan 6 SW ± iR~ 3 ' u [l%Io- fi , 3 V co W 0 \oj\ l/i . (42 b) 

COS / v sw 
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Because T^o is a complex constant, the viscous effects not only modify the inviscid growth 
rate behaviour but also introduce an 0(R~ 3 / 14 ) correction to the inviscid result indicated by the 
leading term in (42b). As ft 0 — * oc, the viscous effects become weaker relative to the inviscid 
contribution, and we recover the earlier result (27). On the other hand, as fto — > 0, the viscous 
term dominates the eigenvalue correction to yield 

« ~ -R-Wfi otanfl™ ± ifi- 3 / 14 ^ 73 S1 ^7 [-KoW 0 '(0)] 1/2 . (43) 

With the proper choice of sign for the second term on its right-hand side, eq. (43) yields a 
vortex mode that amplifies in the X direction, but at a growth rate that is relatively small 
in comparison with the real part of a. However, because the rate of decrease in the growth 
rate as fto — ► 0 is significantly slower than that of the real part of a, the eigenvalue-correction 
term in (43) will eventually become as large as the leading term. This indicates that a different 
asymptotic structure must be considered at this stage. Not surprisingly, the onset of this new 
asymptotic structure corresponds to fto = 0(/?~ 3 / 28 ) or, equivalently, ft = R~ x ^ 4 ftoo [ftoo = 0(1)], 
which is the triple-deck scaling. 

4.2 Viscous-inviscid Interactive Regime 

For ft — R~ x / 4 ft oo, the disturbance structure is fully viscous-inviscid interactive in nature, and 
the perturbation quantities within the lower deck expand in the form 

(■ u,v,w,p) = (Uo,R- 1,2 Vo,W 0 ,R-' /4 Po) + - - ( 44 ) 


This expansion leads to the following set of governing equations after substitution into the 
continuity and momentum equations; 


*(«ooCo + AhM) + ^ = 0, 


dYn 


00 


^ [ooo^os 0 SVJ u'(0) + P oosin 0 sw u/(O)] >oo Co + cos 0 su ,u'(O) V 0 

dP , b 


dYr 


0, 


00 


i [«oocos 6 sw u'( 0) + /?oosin 0 su ,u/(O)] >oo W 0 + sin 0 4U/ u/(O) Vo = 


• P+ d ' 2Uo 

-i otoor o + 2 » 

Cl II 


, „ „ (P w 0 

-IP 00P0+ -JyT- 


(45a) 

(45 b) 
(45c) 
(45 d) 
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Here, 

is the new lower deck coordinate; 


y 00 = rr'"i 


(45e) 


Poo = R' h P, ( 45 /) 

is the scaled (i.e., 0(1)) wave-number parameter; the corresponding, scaled eigenvalue tt 0 o is 
defined as 

«oo = R y4 <y- (45 g) 

The boundary conditions associated with (45a-d) are given by the no-slip conditions 

[J 0 = Vq = Wq = 0 at the surface Vqo = 0, (46a, b, c) 


with the requirements 

Uq ~ cos 0 su ,a'(O) A, Wo = sin 9 sw w'(0) A as — >■ oo 


(46c/, e) 


for matching with the inner limits of the main-deck solutions. The displacement-thickness per- 
turbation A in (46d, e) is related to the pressure perturbation Pq via the interactive relationship 


(a 00 cos 9 

SUf + dposin 9 SW ) a 

i 0 — ~ A" 


(46/) 


(«oo + dlof 2 

The triple-deck problem defined by (45a-d) and (46a-f) is equivalent to the zero-frequency limit 
of the problem considered by Smith and Stewart (1987) who analyzed the effect of crossflow on 
the behaviour of Tollmien-Schlichting (TS) instability. Accordingly, the dispersion relation for 
the interactive stationary vortices can be easily recovered by setting the frequency parameter 
equal to zero in equation (Bl) of Smith and Stewart. This relation is given by 


(«w + Poof 12 (ttoocos 0 SW + /? 00 sin 9 sw f + 3 Ai'(O) {/[a 0 ocos 0 su ,u'(O) + /?oosin 0 S uV(O)]} :>/ ' , = 0 . 

(47a) 

The dependence of the normal-mode properties on the shape of the base velocity profile has 
disappeared completely because the leading order eigenvalue cv 0 o can be determined with a 
knowledge of just the wall-shear components u'(0) and u/(0). 

Two of the roots of the dispersion relation (44) behave in accordance with 

, 1/3 sin0 su , f T/ ,, ml i /2 


5/3 _ 


c*oo ~ — Poo tan 9 SW i 


oo 


cos 1 / 2 9 S 


— KooW’o^O)] ' when j3 00 oc . 


(47b) 
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Hence, these match onto the stable and unstable modes from equation (43) when 1 << (3 00 << 
^> 3/28 Q n t j ie Q ti ier hand, numerical calculations (not shown in this paper) have revealed that 
as ft oo — ► 0 the analytical continuation of the unstable mode in (47b) behaves in the manner 


a . „ 9 / 5 |cos0 su , - tan0 u ,sin0 su ,| 6/5 

a m ~ /9oo tan 0 W tfto A [-:Mi'(0)] 3 / 5 cos 8 / 5 


(47c) 


where 0 W = arctan[tan 0 su ,|t^] and A = [cos 2 # stu u'(0) 2 + sin^^w^O) 2 ] 1 / 2 define the direction 
and the magnitude, respectively, of the wall shear stress vector [cos sin ^^'(O)]. Equa- 

tion (47c) implies that the mode that is unstable at large values of ftoo also remains unstable at 
very small values of ftoo, although the spatial growth rate decreases at a faster than linear rate 
when ftoo — ► 0. We also made an attempt to solve the dispersion relation (47a) numerically in 
order to track the behaviour of the stable mode in (47b) across the ft 0 o = 0(1) range and into 
the ftoo « 1 regime; however, in each case, the Newton- Raphson method failed to converge 
after ftoo became sufficiently small. Consequently, the behaviour of the stable mode could not 
be determined either analytically or numerically. 

ft is worth noting that the large wave-number limit implied by (47b) also represents nearly 
neutral vortex modes (as in the small wave-number limit (47c)). However, the large wave-number 
modes are aligned with the local free-stream direction, whereas the small wave-number modes of 
(47c.) are aligned with the local wall-shear direction. Substitution of (47c) into the momentum 
equations (45b, d) implies that the convection term associated with the linear behaviour of the 
local directional mean profile tends to zero as ftoo — 1 ► 0. Of course, this decrease in the magnitude 
of the convection term cannot continue indefinitely because when ftoo = 0(R ~ 5 / 28 ) ) or 


09 = 0(1)), 


(48a) 


one must consider the more refined approximation for the mean-flow profiles 

u ~ U\Y + u{Y 2 + ... (as Y — * 0) 


(486) 


and 

w ~ w^Y + W 2 Y 2 -f ... (as V — ► 0) (48c) 

rather than the linear approximation that has been sufficient thus far. For brevity of notation, we 
have denoted the wall-shear components u'(0) and u/(0) by u i and ufi, respectively, in equations 
(48b, c). 
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The solutions for the lower deck perturbations in (44) suggest that as /5 0 o — > 0 the component 
of the horizontal velocity perturbation along the direction of the wave-number vector (i.e., 
the component that is proportional to the divergence of the horizontal velocity field (f/ 0) W 0 )) 
becomes asymptotically small in comparison with the component that is orthogonal to the 
wave-number vector (and, hence, is proportional to the vertical vorticity field). However, the 
dispersion relation for the unstable crossflow mode is still determined by the former component 
which, along with the perturbations in vertical velocity and pressure, is decoupled from the 
vertical vorticity field. In view of this decoupling, we limit our attention to the former set 
of quantities (in particular, Vo and Po) for ft « R~ 1//4 . The vertical vorticity field and the 
dominant velocity perturbations can be evaluated on an a posteriori basis, if necessary. 


4*3 Interactive Wall- Shear- Aligned Vortices 

For ft — <9(i?“ 3//7 ), a modified interactive structure is obtained for which the eigenvalue a and 
the lower deck perturbations (f/, V, IT, P ) expand as 

a ~ — R~ 3 ^ 7 ft tan 6 W + RT a ^ oq (49a) 


and 

(u, V , w, p) = (Co, R- S/ 7 V 0 , Wo, R~ 3 / 7 Po) + ..., 

respectively; the new lower deck coordinate is 

Y = R 1 / 7 Y. 


(496) 


(49c) 


The momentum equation along Y reduces to a form that is analogous to (37c) and (45c). After 
a suitable manipulation of the continuity equation and the momentum equations along the X 
and Z directions, one finds that Vo is governed by the third-order equation 


d3y ° i(XY + Af 2 )^ + *(A + 2\Y)V 0 = /3 2 sec 2 0 W P O 


dY 3 v ' dY 

and is subject to the wall boundary conditions 


and the condition 


dVp 

dY 


r 0 (o) = ^(h = 0 > = o 

v ’ dY 


-i(A + 2Ay)^ as Y^oo 


(50) 


(51a, 6) 


(51c) 
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for matching of V 0 with the main deck solution 


= -iR - 3 ' 7 


(5k/) 


Here, 


A = aicos0 sw ui, A = //(sin 0 sw w 2 — tan 0 W cos^ su ,m 2 ) (51e,/) 

and the displacement thickness perturbation A is related to the pressure perturbation Pq via 
the interactive relation 


P Q = ftcos 0 w (sin 9 SW - tan 6 W cos 6 aw f A, (51g) 

which is obtained from (46f) and (49a). Note that the set of equations (50) through (51c) can 
be simplified through an affine transformation so that the independent parameter A appears 
only in (51c). However, the advantage gained through this procedure is relatively small; hence, 
we will not introduce this transformation in the present analysis. 

For a given value of ft, solutions to the inhomogeneous boundary value problem (50) and 
(51a-c) are possible only for a limited set of values for A, each of which determines the growth 
or decay rate of the respective eigenmode in view of the definition (51e) of A. In the limit 
of ft — > oo, Aj increases at a faster than linear rate (see (45b)). Therefore, the quadratic 
terms in (50) become relatively small, thereby yielding the required overlap with the governing 
equations (37) and (38) in the triple-deck regime. On the other hand, when ft is reduced, the 
growth rate (as well as the real part of «i) is expected to decrease monotonically; therefore, the 
linear convection term in (50) that involves A will diminish in importance in comparison with 
the quadratic part. As a result, the eigenvalue parameter A will drop out of the leading-order 
equation for V Q . The latter does not contain a pressure-gradient term and is identically satisfied 
by the small-F behaviour of the main deck solution (51d): 

v ~ -iR~ 3/7 AY 2 A, (52) 

which also satisfies the boundary conditions (51a-c) above. The correction to this inviscid 
solution is driven by both the residual convective acceleration that involves A and the inviscid 
solution (52), and also by the pressure-gradient term in the original governing equation (50). A 
balance of the two forcing terms provides the asymptotic scaling A = 0(ft 5 ^ 2 ) for the eigenvalue 
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parameter A in the limit /3 — ► 0. Accordingly, we introduce the rescaled eigenvalue 


1 = /H/ 2 a, 


(53a) 


the rescaled wall-layer coordinate 

? = p' / 4 Y, 

and the scaled pressure-gradient parameter 

A = /3 -1 A - sin 0 sw w 2 — tan 9 W cos 0 SW u- 2 . 
We expand Vo and Po in the forms 

y 0 = -p' /2 i~\Y A + /3 9 / 4 Voi + ... 


(535) 


(53c) 


(53d) 


and 

Po = /3Poi + ... (53e) 

as (i — > 0, where Voi is now a function of Y. If equations (53a-e) are substituted into (50) 
through (51c), Voi is governed by the equation 


d 3 V 0 1 stfdVr 


^01 
FT 
dY 


+ 2tAf Vox = fih \f(h = sec 2 0 w P m + Ih Y A 
dY L 


(54a) 


subject to the boundary conditions 


MO) = = 0 


and 


and the interactive relation 


dV, 


01 


dY 


dY 


— iXA as y — > oo, 


(546, c) 


(54 d) 


P 01 = cos # w (sin 6 SW — tan 0 W cos 6 SW ) 2 A. 


(54c) 


Solutions to this inhomogeneous boundary value problem are possible only if the eigenvalue 
correction A satisfies the solvability condition 


r 

Jo 


fF dY = 0, 


(55) 
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where F represents the eigensolution of the adjoint problem 

- id? Vo, = 0; F(0) = 0, ^£,^£^0 as f oo. (56 a,b,c) 

d dY dY dY 

We can sliow that F is given by 

F = sU(0,s) (57 a) 

within a constant multiplicative coefficient; here, 

/ :\l/4 ; 

5 = UiAj Y (576) 

(with the assumption that the principal branch has been chosen for the multivalued function 
on the right-hand side) and (7(0, s) denotes the parabolic cylinder function (Abramowitz and 
Stegun, 1964). Based on the property that / 0 °° s 3 U(0,s ) = (7(0,0), the dispersion relation in the 
0 « 1 regime is given by 

\ n fl t n n V 2 [/o°° -8(7(0, s)rfsl (iA) 1/2 , 7 ) 

X — -sec. 0 w (sin 9 SW - tan 0 W cos V sw ) jjj — ^ j ■ \p(c) 

The numerical value of the bracketed ratio (which involves parabolic cylinder functions) on 
the right-hand side of (57c) was quoted by Hall (1986) to be 1.198, approximately, and our 
calculations verified this result. In view of (49a), (51e), and (53a), equation (57c) implies a 
positive growth rate for the crossflow vortex under consideration. Equation (57c) is, in fact, 
identical to the dispersion relation obtained by Hall (1986) in the context of a rotating-disk 
flow, provided that the Coriolis terms are set equal to zero in the latter equation. 

Equations (53a) and (57c) imply that for R~ T « 0 « 1 (r > 0) the spatial growth rate of 
the crossflow vortex, as well as the wave-number correction given by the real part of , decreases 
at a rate that is proportional to 0 5/2 . A comparison of (57c) with (48b) reveals that this rate 
of decrease is significantly greater than the corresponding decrease (shown in section 4.2) when 
0 was reduced below the range 0 = O)/? -1 ^ 4 ), which corresponds to the triple-deck regime. 
This decrease in is accompanied by a corresponding decrease in the wave-number-aligned 
component of the disturbance slip velocity at the inner edge of the main deck. (See (54d).) 
Note in (51a) through (54e) that this latter component gives rise (through viscous retardation) 
to the displacement perturbation A within the lower deck. The displacement perturbation 
induces a pressure disturbance in the upper deck that is fed back to the lower deck, which 
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sustains the interactive structure of the disturbance. As a consequence of the decrease in this 
slip-velocity component as ft decreases, the higher order corrections (which result from the 
previously neglected cubic and higher order derivatives of the directional mean profile near the 
surface) exert an increasing influence on the outer behaviour of the correction term Voi within 
the lower deck. For instance, when ft = 0(/? _1/14 ) (i.e., ft = 0(f?~ 1/r2 )), the boundary condition 
(54d) must include the effect of the cubic term in the small-T Taylor expansion of the main-deck 
solution (51d). (See Hall, 1986.) Unlike the inviscid leading-order solution in (53d) for which 
the viscous-diffusion term is equal to zero, the inviscid contribution associated with the cubic, 
term does not satisfy the viscous lower deck equation. The viscous correction to this inviscid 
contribution would modify the viscous displacement function A by an 0(1) amount and hence, 
also, the dispersion relation (57c.) that determines the growth rate of the vortex. However, by 
using the mean boundary-layer equations, the third-order derivatives of u and w can be shown to 
be equal to zero for a large class of boundary-layer flows, such as boundary layers with no suction 
at the surface. In that case, the dispersion relation (57c) will remain valid until ft = 0(R~ 8 / 63 ) 
(i.e., ft — 0(/?~ 5 / 9 )), at which stage the fourth-order derivative of the directional mean profile 
(= 24(sin O sw w 4 — tan 0 W cos O sw u 4 )) will come into play. 

The cubic (or quartic) terms can reduce the growth rate of the crossflow mode to negative 
values in the ft — 0(/? _1/2 ) (or ft = 0(R~ 5 ^ 9 )) regime and cut off the crossflow instability at 
sufficiently small wave numbers. However, at least for the rotating-disk case Hall’s analysis 
(Hall, 1986) indicates that these terms do not significantly affect the stability properties. Al- 
though these aspects are interesting and important in regard to their influence on the receptivity 
properties (section 5), we will not address them in this paper. 

5 Concluding Remarks 

The stability properties of long- wavelength stationary crossflow modes were examined, beginning 
with the long-wavelength Rayleigh regime R~ T « ft « 1 (r > 0). In the latter limit, the 
leading-order solution to Rayleigh’s equation reduces to a simple, vertical displacement of the 
jetlike crossflow profile. A nominally 0(ft 2 ) correction term (that involves the operator Ai 
in Rayleigh’s equation (4a)) gives rise to a significantly larger amplification rate, specifically 
0(/?U 2 ) relative to the real part of the wave number o. This effect is caused by a commensurate 
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magnification across the boundary layer in the magnitude of the singular homogeneous solution 
(viz., vb ) to the long-wavelength form of Rayleigh’s equation. The expression for the leading- 
order growth rate resembles to an analogous expression derived by Drazin and Howard (1962) 
for the sinuous mode of a symmetric, jet or wake profile. The reason for this similarity is that in 
both cases the above-mentioned displacement solution satisfies both of the imposed boundary 
conditions. The most striking difference between the long-wavelength stability properties of the 
boundary-layer flows considered in this paper and the free-shear flows considered by Drazin and 
Howard (1962) results from the earlier onset of viscous effects in the boundary-layer case. In 
particular, as /? is reduced below the long-wavelength Rayleigh regime, the viscous effects lead 
to a cascade of asymptotic subregimes across which the orientation of the crossflow vortices 
changes from the direction of the local inviscid streamline to that of the mean wall-shear vector. 
The variation in the spatial amplification rate as /? varies through the subregimes is shown 
schematically in Figure 4. 

The study of the effect of surface curvature in section 3 demonstrates that the convexity of 
the surface along the longitudinal direction can stabilize the inviscid crossflow modes. A similar 
effect can be expected in the viscous regimes as well. The influence of convex surface curvature 
is particularly prominent near the nose of an airfoil and, indeed, has been noted in a number of 
investigations including those by Malik and Poll (1985) and Collier and Malik (1987). 

We now close the discussion by commenting on other possible extensions and applications of 
the present theory. One important extension of these results is the consideration of compressible 
base flows and nonzero disturbance frequencies. Both generalizations require relatively minor 
modifications to the theoretical framework of sections 2 through 4; however, the modifications 
promise results of considerable practical importance. In particular, the inclusion of a nonzero 
frequency parameter would indicate how a lower branch of the neutral stability curve emerges 
in the wave-number regime studied herein as the frequency parameter increases to the Tollmien- 
Schlichting (TS) regime considered by Stewart and Smith (1988). The development of this 
analytical connection between crossflow and the TS modes as the Mach number is increased 
to supersonic values will be interesting. The TS instability of three-dimensional supersonic 
boundary layers has previously been considered by Choudhari and Streett (1990) and in further 
detail by Ryzhov and Terent’ev (1991). Other relevant works on unsteady crossflow vortices 
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include the paper by Bassom and Gajjar (1988) on high-frequency modes and the work of 
Balakumar, Hall, and Malik (1991) on the algebraic growth of low-frequency modes in a rotating- 
disk boundary layer. A comparison of the results of the abovementioned extension with the 
results obtained by these authors will be useful. The effects of lateral and longitudinal curvatures 
in the supersonic regime may also be studied by combining the present analysis with that of 
Duck and Hall (1989) and Dando (1992), respectively. 

Another extension of the present study that we are currently pursuing concerns the wave- 
number regime of j3 < 0(R~ 1 / 2 ), where a lower branch of the neutral stability curve may 
exist for stationary- vortex instabilities. The possible existence of a lower branch of the neutral 
stability curve for stationary crossflow vortices on a relatively flat, nonrotating surface has im- 
portant practical ramifications, in particular for receptivity (Choudhari, 1994). The absence of 
a lower branch would imply that the induction of crossflow vortices near the leading-edge region 
(whether by surface disturbances or by free-stream modulations) is likely to be the most impor- 
tant source, overall, of stationary vortex modes in a swept-wing boundary layer. Secondly, the 
nature of the receptivity caused by nonlocalized, stationary surface roughness will be consider- 
ably different from a similar mechanism in the case of nonstationary modes (Choudhari, 1993). 
Choudhari (1993) showed that the receptivity due to an elongated region of small-amplitude 
surface nonuniformities can be evaluated as an integral over each infinitesimal subregion. In 
each of these subregions, a suitable extension of Goldstein (1985) and Ruban’s (1985) theory 
can be applied to predict the local receptivity. When the neutral stability curve possesses a 
lower branch, the behavior of this integral is similar to that of a standard Laplace integral if 
the distribution of surface nonuniformities is periodic. The saddle point of the above integral 
coincides with an interior location that corresponds to a near resonance between the wave num- 
bers of the periodic forcing and the instability mode (Zavol’skii et al. 1983). Thus, the most 
significant receptivity is produced by a narrow range of locations in the vicinity of this saddle 
point. In the absence of a lower neutral branch, on the other hand, the above integral becomes 
dominated by end-point contributions; this indicates that the receptivity is dominated by the 
far-upstream nonuniformities. Thus, the dependence of receptivity on the roughness geometry 
in this case must be compared with the simple behaviour observed by Zavol’skii et al. (1983) and 
Choudhari (1993, 1994). Finally, we note that the efficiency of receptivity caused by a localized 
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roughness element in the ft « R~ 1 ^ 4 regime is also expected to be quite different from the 
analogous localized receptivity problem in the case of the TS waves studied by Goldstein (1985) 
and Ruban (1985). This difference is because a small-amplitude roughness element produces 
very little disturbance in the vicinity of the direction that is orthogonal to the wall-shear-stress 
vector and the latter direction coincides with the direction of the wave-number vector for the 
crossflow-vortex modes in the above wave-number regime. 

In section 4, we indicated how the dispersion relation from one asymptotic regime matches 
onto the appropriate limit of the dispersion relation from an adjacent region when the relevant 
wave-number parameter is allowed to either decrease to zero or tend toward infinity. Physically, 
this matching procedure is equivalent to considering a fixed location on the body surface and 
varying the wavelength of the crossflow disturbance. This procedure was done partly for want 
of a convenient physical model for a three-dimensional boundary layer; however, in practice, a 
more relevant approach may be to explore the results of an analogous matching procedure in 
which the disturbance wavelength is kept fixed but the Reynolds number is varied by moving 
toward the nose region along the body surface (so that the mean boundary-layer profiles must 
also vary accordingly). Because the long-wavelength modes are nearly neutral (except within the 
triple-deck regime), it will be useful to examine whether any modulations of the incoming flow 
along the z direction can directly excite the unstable vortex modes near the nose of the body. 
Problems of this type have not yet been studied in the context of three-dimensional boundary 
layers, although some interesting results have been reported for the case of a two-dimensional 
boundary layer by Goldstein, Leib, and Cowley (1992) and Goldstein and Leib (1993a, b). A 
similar receptivity problem that involves the generation of Gortler vortices has been studied by 
Hall (1990). 

Great strides have been made recently in the analytical description of the nonlinear stage 
during the transition process. (See, for example, the milestone works by Hall and Smith (1991), 
Goldstein and Leib (1988), and Goldstein and Choi (1990).) However, with the exception of a 
recent effort by Bassom and Hall (1991b), not much research has been devoted to the nonlinear 
development of crossflow-vortex instabilities. Indeed, because the initial amplitudes of stationary 
crossflow vortices can at times be much larger than those of nonstationary modes (Choudhari 
and Streett, 1990), the nonlinearity is likely to significantly influence the development of the 
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stationary vortices much earlier than in the case of the nonstationary vortices. Because the 
amplification rates for the stationary vortices remain relatively small in a significant part of the 
long-wavelength regime, the nonlinear effects might have an influence within the same regime. 
The likelihood of this possibility in flows of practical interest must be determined. 
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Figure 1 . Schematic of asymptotic structure of long-wavelength Rayleigh modes. 
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Figure (2b). Long-wave length asymptotes of curves in figure 2(a) 
(-Real(a) and -10*Imag(ot) as functions of wave-number 
parameter S for (3 H = 0.6 and Q — 45°. 







Figure 2(d). 
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Comparison of long-wavelength asymptote (34) 
with numerically computed neutral curve 
(Blackaby and Choudhari, 1993) for stationary 
vortex disturbances in three-dimensional 
boundary layer over convex surface. (Upper 
branch asymptote of Gregory, Stuart, and 
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Figure 4. Asymptotic subregimes in stationary crossflow instability in three-dimensional boundary layer over flat surface 
(sketch not to scale) where a, denotes imaginary part of the eigenvalue parameter a, and $ mi , and represent wave 
numbers that correspond to maximum spatial growth rate and upper branch of neutral stability curve, respectively. 
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